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Abstract. We investigate linear and nonlinear transport in interacting single- wall 
carbon nanotubes (SWCNTs) that are weakly attached to ferromagnetic leads. For 
the reduced density matrix of a SWCNT quantum dot, equations of motion which 
account for an arbitrarily vectored magnetisation of the contacts are derived. We 
focus on the case of large diameter nanotubes where exchange effects emerging from 
short-ranged processes can be excluded and the four-electron periodicity at low bias 
can be observed. This yields in principle four distinct resonant tunnelling regimes, but 
due to symmetries in the involved groundstates, each two possess a mirror-symmetry. 
With a non-coUinear configuration, we recover at the 4N 4N ± 1 resonances the 
analytical results known for the angular dependence of the conductance of a single 
level quantum dot or a metallic island. The two other cases are treated numerically 
and show on the first glance similar, yet not analytically describable dependences. In 
the nonlinear regime, negative differential conductance features occur for non-coUinear 
lead magnetisations. 



PACS numbers: 73.63.Fg, 85.75.-d, 73.23.Hk 
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1. Introduction 

Carbon nanotubes are promising candidates for constituent of tomorrow's electronic 
nanodevices [1] and thus subject of today's theoretical [2], [3], HI O [6l [3, [8] and experi- 
mental investigations P[10l[IIlll2l[l3l[llll5lll6l[l7l[l8]. A single- wall carbon nanotube 
(SWCNT) is just formed by the planar honeycomb lattice of graphene, wrapped along 
one axis such that a closed, seamless cylindrical surface arises [19]. To guarantee this, 
the axis and the tube diameter cannot be chosen completely arbitrary, but in princi- 
ple there are countless possible configurations. Only special ones, however, inherit the 
metallic properties of graphene, which stem from the fact that in graphene conduction 
and valence band touch at the corner points of the first Brillouin zone. Two of those 
points are independent and labelled Fermi points ±Kq. To be metallic, a SWCNT must 
contain the Fermi points in its reciprocal lattice of allowed momenta and all armchair 
type tubes, as well as certain types of zig-zag configurations [19] fulfil this condition (but 
curvature effects yield a small band gap for all SWCNT types except for the armchair 
ones). For low energy processes, just momenta close to the Fermi points can contribute, 
and as the SWCNT diameter is normally a small fraction of the tube length, the size 
quantisation permits many longitudinal, but merely one radial mode. This fact makes 
the carbon nanotube a one-dimensional conductor. Moreover, the shape of the graphene 
bandstructure is very beneficial in the considered region: it develops linearly around the 
Fermi points, which allows to describe the electronic properties of a nanotube by means 
of the Tomonaga-Luttinger model for interacting fermions [21 [5], involving a bosonisa- 
tion of the electron operators. 

A research topic in rapid expansion is spin-dependent transport. In particular, recent 
theoretical works have focussed on transport across interacting quantum dots [2011211 [22] . 
metallic islands [23] and wires [5] with non-coUinearly magnetised leads, offering the 
prospect of spin-sensitive single electron transistors. It has been already realised that 
interactions [2011211 [23] as well as reflection processes at the lead-system- interface [2^ [22] 
can strongly influence the spin-accumulation in the dot (island), and hence the properties 
of the spin-valve transistor. During the last years, there have been various experimen- 
tal investigations on spin- dependent transport in SWCNTs [HI [13], [TU [T71 [18], and a 
possible measurement setup is sketched in Fig. [H Notice that we include an arbitrary 
relative magnetisation of the leads and do not restrict ourselves to a common limitation 
of nowadays experiments on SWCNTs, where merely the two discrete states of parallel 
(P) and antiparallel (AP) contact magnetisation can be realised (0). There have already 
been experimental studies on non-coUinear spin transport in magnetic multilayers [25] . 
so we can expect that in the near future it will be possible to attach arbitrarily polarised 
contacts to a SWCNTs. 

I To be able to obtain an AP configuration, usually source and drain contact are manufactured with 
different width so that the coercive fields emerging in the metals are not of equal strength. Performing 
a sweep of an external magnetic field, the wider contact switches later to the opposite magnetisation 
than the smaller one and for some range of the applied field the leads are thus polarised in antiparallel. 
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Figure 1: Ferromagnetic contacts are attached to the SWCNT. Besides a bias voltage Vbias, a 
gate voltage Vgate can be applied via a back gate in order to shift the chemical potential inside 
the tube. The magnetisations mg and rhd are controlled by varying an external magnetic field. 

A variable of special interest in spin-dependent transport experiments is the 
so-called tunnelling magnetoresistance (TMR), which we define for coUinear contact 
magnetisations as (TMR) = {Gp — Gap)/Gp, where Gi is meant to denote the 
conductance in the parallel {i = P) resp. antiparallel {i = AP) configuration (§|). 
In Ref. [I7j, the observation of quite regular TMR oscillations with the gate voltage is 
reported, where large changes of almost 20% in the conductance are reached. Moreover, 
the TMR acquires negative values, which means that the AP conductance can exceed 
the P one. 

So far, spin-polarised transport in interacting SWCNT quantum dots has been 
considered only for very long tubes, which do not exhibit level quantisation and 
charging effects [5]. For this latter case, a non-trivial dependence of the current on 
the interaction strength was predicted. Recently, a single impurity Anderson model 
with four degenerate orbitals was proposed [22] as a minimal model to understand the 
aforesaid negative TMR effect reported in [17] . Particularly, the authors could show that 
it can originate from multiple reflection processes at the SWCNT-contact interfaces. 

In this work, a microscopic treatment of spin-dependent transport in SWCNT quantum 
dots is presented. Specifically, we focus on the low transparency regime, where a weak 
coupling between tube and contacts is assumed. Therefore it is justified to treat all 
Hamiltonians associated with the tunnelling barrier in lowest nonvanishing order. The 
opposite case of high transparency, i.e. low ohmic contacts, has already been studied 
both experimentally [18] and theoretically [26l [27]; it was found that phase shifts 
picked up during backscattering events at the tube ends yield a Fresnel/ Fabry- Perot 
interference pattern in V^ate — Hias— TMR plots . 

§ There are various other definitions of the TMR, e.g. {Rap - Rp)/Rp [H HI], (Gp - Gap)/{Gp + 
Gap) [22]! '^{Rap — Rp)/{Rp + Rap) [H], where Ri denotes the resistances in the parallel {i ~ P) 
resp. antiparallel {i = AP) configuration. The physical meaning, nevertheless, is always the same, as 
merely the normalisations differ. 
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Our model takes into account interface reflections, as well as virtual transition processes. 
Both are relevant exchange effects for spin-dependent transport |23], inducing a 
precession of the spin accumulating on the quantum dot. 

For strict lowest order perturbation theory, nevertheless, the interface reflections cannot 
be source of any negative TMR, as it will become clear in the course of this article, 
which is structured as follows: 

In section [2] we describe the model Hamiltonian we use for our system, section [3] 
explains how to derive the master equation for the SWCNT density matrix in the 
presence of arbitrary lead magnetisations. Section |4] introduces the necessary coordinate 
transformations to proper spin quantisation axes and subsequently, section [5] contains 
the results we acquire for the current in the linear transport regime. Section [6] presents 
numerical data for the nonlinear case and finally section [8] gives a summary of our 
achievements. In the appendix, some explicit calculations skipped in section [3] can be 
found. 

2. Model Hamiltonian 

To build up a model for spin-dependent transport in the quantum dot regime of a 
SWCNT, we expand the theoretical work [8] on correlated transport in carbon nan- 
otube quantum dots to include spin polarisation of both contacts. The magnetisations 
rhs and fhd may enclose an arbitrary angle 6, which is a possibility future experiments 
are likely to offer (see Fig. [1]). 

The Hamilton operator we use for the setup in Fig. [1] is 

H = Hq+J2Hi + Hti + Hm , (1) 

l=s,d 

where effects of external voltages have been absorbed into Hq and Hi. The index / labels 
the source {l=s) and the drain {l=d) contacts, which are metallic and thus characterised 
by a Hamiltonian 

= 5Z ( ^9 - Sgn((Ti)Estoner " cVJ) cf^^ .Q,,,^- . (2) 
Q 

In (|2]), Qo-jg- is the fermionic annihilation operator for an electron of momentum q in 
lead I with spin a/ G {+/, — where ±; denotes the majority /minority spin species for 
a quantisation along ai. For reasons of simplicity, we just apply the Stoner model for 
our ferromagnetic contacts, but any other description could be used as well. Both the 
bias voltage VJ and the Stoner exchange splitting energy -Estoner for the two spin species 
in the ferromagnetic leads add a shift to the kinetic energy e^- of the particle, e is the 
electron charge. 
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The tunnelling processes at contact / are modelled by 

Hti = J2 [ (^'(^l ^Ui^^i^^i^ + h.c.) . (3) 

The so-called transparency Tiif) specifies the tunnelling properties of lead /, '^QfJ^{r) 
and \E'/o-((?^ are the electron annihilation operators in real space representation, for a 
particle of spin cx; in the tube resp. in lead I. For later purposes it is necessary to know 
the decomposition 

ili^Xr) = Jde Di^XE\:il)Y.^i<i{^ (4) 

where DlfJ^ is the density of states for carriers of spin o"; in lead /, taking as an argument 
the total energy 

^tot If = e-eVi- Sgn(cri)£^Stoner 

of a particle, (piqir) is just a wave function inside the lead. 

Furthermore, we allow for boundary backscattering by introducing the 'refiection' 
Hamiltonian 

Hm = - I dV AKr) 5^sgn(aO^U(O^0.,(^1 , (5) 

which is equivalent to the momentum space expression used e.g. in [23] (where also 
the relation to the mixing conductance introduced in [28] is explained). An electron 
picks up some phase when it is scattered at the tube ends, which overall results in a 
certain energy shift, coming with a positive or a negative sign depending on the spin 
polarisation. In other words, being close to a contact, the electron feels a bit of the 
magnetic field inside it, causing a position dependent energy splitting A/(r) for the two 
spin species. 

We assume the tunnelling barrier to be spin-independent and as there is only weak spin- 
orbit coupling in clean nanotubes, it is justified not to consider any spin-fiip processes. 
At this stage, we still want to stay general and do not introduce any assumption on the 
position dependence of Ti{r) or Ai{f). Later on, however, we will have to impose the 
restriction that both parameters are of relevant value only nearby the contacts. 

The most complex contribution to the system Hamiltonian is given by the terms which 
belong to the SWCNT itself, due to the presence of Coulomb interaction. 
The SWCNT Hamiltonian Hq = Tq + Vq reads 

Hq = hvp ^ sgn(f ) ^ n ct^g„ Cf^^n + 

+^ E // d\'¥^,^{T)¥^^,(f')V{T- r')^0.^(r')^0.e(^1 • (6) 
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Notice that we index the electron operators by a spin cr0, cTq G {Tj i} which refers to an 
arbitrary, but fixed unique quantisation axis inside the tube as e.g. sketched in Fig. |H 

The first term in ([6]), Tq, collects the kinetic energy of particles in the nanotube. To 
recognise this, one needs some knowledge about the SWCNT bandstructure. We already 
pointed out in the introduction that from the graphene lattice, the reciprocal space of 
an armchair SWCNT inherits two Fermi points ±Ko, where valence and conduction 
band touch. Transport can take place in the vicinity of these points, around which the 
bands develop linearly and, due to our restriction to low energies, one-dimensional. Fig. 
[2] (left) shows the two branches we obtain at each Fermi point when periodic boundary 
conditions (PBCs) are employed [I9] : a left mover band L of negative slope and a right 
mover band R of positive slope; both slopes have an absolute value hvp, where vp is 
the Fermi velocity. By the linear transformations 

VL[~K„][^Ko]{r)) , 

(7) 

wave functions (pfKn which fulfil open boundary conditions (OBCs, Fig. [2] (right)) - those 
are the appropriate ones for the finite-size system the tube represents- are constructed 
from the usual PBC wave functions ^tk^f ij ^ {L, R} )■ The latter can, as worked out 
in [8], be freed from the dependence on the momentum Kn by the approximation 

^rn^F{r) ^ e'"""</^.F(r) . (8) 

The transformation ([7]) maps the four rF branches onto two bands L and i?, whereas 
the number of admitted momenta Kn = hv fU , n G Z doubles [8]. 




periodic boundary conditions (PBCs) 




open boundary conditions (OBCs) 



Figure 2: Band structure of a SWCNT with periodic (left, PBC) and open boundary 
conditions (right, OBC). For (non-interacting) wave functions which fulfil OBCs, the spectrum 
is characterised by two possible bands L and R and momenta Kn = hvfn , n ^7L. The energy 
levels of left- and right mover branches can be shifted with respect to each other, which results 
in a nonzero band offset b. 
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The second part of Vq, arises from the strong electron-electron interactions inside 
the nanotube, where V{f—f') is the Coulomb potential, including screening. Under the 
assumption that forward scattering processes are the only relevant ones, this potential 
part can be simplified enormously. We are allowed to do so if interactions within one of 
the two graphene sublattices, so-called intra-lattice interactions, cannot be distinguished 
from inter-lattice interactions. Such an approximation holds for tubes with sufficiently 
large diameters ( > 5 lattice constants) and length ( > 100 nm) |29j . 

By combining the OBCs with the expression ([8]) and the linear character of the bands 
around the Fermi points, the electron operator 

f£{L,R} " 

can be expressed in terms of the ID operator i^fa-Qpix) = -^7= e'^^'^^'^^'^"^ CfaQn as: 
^Q.Jf) = v^^sgn(F)5, 

rfF 

with r G {L, R}, +r = r, -r = f {L = R, R = L), 6r,f ^ if (r, r) G {(L, L), [R, R)} 
and Lt being the SWCNT length. 

Introducing PfaQpix) := "4^1^ . Fi^)''i'faQFix) , Vq can be rewritten as: 



ro-Q 



Ve 2 



1 V r^t r^t 

2X1X1X1 / / (ixdx'pfaQF{x)V{x,x')pf'a'^^F'{x') . (10) 



ff FF a0<7Q 



The one-dimensional density operator PfaQF{x) has the convenient property that from 
its Fourier components, bosonic operators h„^^n , n G can be constructed, where the 
ones indexed with n > stem from ppi^^p, the others with n < from Pi„^p- Within 
the spin-charge-separation model, a unitary transformation to charge- and spin-like 
excitations is usually performed: 

It is possible to reexpress Hq in terms of these new quantities and, by means of some 
more hnear transformations {bjn • • ■ ^ djn , j G {c, s}) eventually diagonalise the 
SWCNT Hamiltonian. 

The final form of = Tq + Vq is: 




eo 



/A/4 - \ 

X ( + ^ Sgn{f)//fae ) + X X ^jria]^^jn ■ (H) 



fa,j band offset n^O j=c,. 



bosonic modes 



The first three contributions flTT]) contains are purely fermionic, as the operator A/fo-Q 
is just defined to count the particles in band fa^ G {L f,L [} and 
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The three summands, as imphed in fill I) , are the charging energy, a term accounting 
for the Pauh principle if more and more electrons are filled in the same band, 
and a correction for a potential band offset S, as the energy levels of L- and 
i?-band might be shifted with respect to each other (illustrated in Fig. [2]). 
Here, eo = nhvf/Lt is the level spacing and Ec = Wqq, where Wnn = 
^ J^^ dx Jq* dx' V{x, x') [cos(?7,{a: + x'}) + cos(n{x — x'})] . For a typical SWCNT, eo 
and Ec are both of the order of some meV. 

The last term counts the energies of collective, bosonic excitations. The energies 
Cm fire for 77. > dependent on the interaction, 6c+|n| ~ |n|eo\/l + 8W„„eo\ and 
thus called charged modes, while the other neutral modes only scale with the level 
spacing: es-|-|„| = ec_|n| = \n\eQ. None of the bosonic excitations influences the particle 
numbers in the single bands and that is why one can classify the eigenstates of the total 
SWCNT Hamiltonian by a vector \N,m) with = [Ni^, N^, Nj^^, Nj^-^) determining 
the fermionic configuration and m counting the bosonic modes. This uniquely fixes a 
state; without any bosonic excitations, no 'holes' are allowed: all Nf„^ electrons in a 
certain band faQ have to populate the Nra-Q lowermost states. 

With this background, we are now ready to start calculating the dynamics of our system. 



3. Equations of motion for the reduced density matrix 

We want to investigate the time evolution of our system consisting of SWCNT quan- 
tum dot and leads by using the Liouville equation for its density matrix p\t) in the 
interaction picture. This representation is well-suited, because we intend to treat both 
the tunnelling Ht = Hts + Hxd and the reflection Hr = Hrs + H^u as a perturbation 
Hi = Ht + Hr to Hq := Hq + Hi. 

Indeed it is a critical question whether or not Hn should be included in Hj. Consider- 
ations in Ref. [23] show that the values of the phase shifts picked up during boundary 
reflections are of the same order as the transmission coefficients, such that a weak con- 
ductive coupling will bring about a weak ferromagnetic coupling. 

In Ref. [22], an Anderson model is used and both the reflection and the tunnelling are 
treated non-perturbatively. 



The equation of motion reads: 



dt 



(12) 



where Hj had to be transformed into the interaction picture according to Hj{t) = 
eiHo{t-to)Hj(to) e-:^^»(*-*") . 

Our flnal interest is dedicated to transport through the SWCNT quantum dot, thus we 
would not mind to lose all the information about the contacts contained in p^{t): it is 
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sufficient to calculate a reduced density matrix (RDM) PQ{t), where the lead degrees of 
freedom have been traced out: 



d 



p^(t) := Trieads {p\t)) , yielding ih-p'^it) = Trie^ds Hj{t),p\t) 



(13) 



3.1. Perturbation theory in the tunnelling and the reflection Hamiltonian 

In general, the contacts can be considered as large systems compared to the SWCNT. 
Besides, in our case of low transparency, the influence of the nanotube on the leads is 
marginal. That is why they can approximately be treated as reservoirs which stay in 
thermal equilibrium all the time, and the following ansatz is valid: 



p\t) = pi{t)p,pd + 0{Ts, Trf) + 0(A„ A, 



(14) 



where pi {I = s,d) are the thermal equilibrium density matrices of the source and the 
drain contact. Putting (HM into (fT3l) . we are allowed to neglect the corrections to the 
factorisation, because in the product with Hj they all produce a higher order remainder. 

Now HJp{t) contains operators belonging to the contacts: either ^!ifj^{r) or ^\^^{r) 
takes part in each term of (see Eq. (|3])), while H^^ (Eq. ([5])) involves no lead 
electron operator at all. Since ps and pd appear under the trace, they give us the 
thermally averaged expectation values of the source and drain operators. Here, due 
to Wick's theorem, all terms with an odd number of lead electron operators, such as 
Trieads [-^t(^) 5 P0(^)P^*Pj vanish. 

Tr leads [-f^ij (i), P0(^)PsPd] is nonzero and of ffist order in the reflection parameters. In 
order to go to the ffist nonvanishing order in the tunnelling, we reformulate (1121) to an 
integral equation: 



p\t) = p'ito) - 7 / dti [^/(ti),p(ti 



h 



to 



(15) 



The time to is chosen to be the time where the interaction is switched on, such that for 
t' < to we have Hj{t') = H^{t') = H^^{t') = and hence p^t') = pQ{t')pspd holds true 
accurately. 

Eq. f|T5l) is exact and we are allowed to reinsert it into f|T2l) . As we are comfortable with 
[Hj^{t), p\t)], we only replace p\t) in the commutator term involving H^{t) and find: 



■T Trieads 

h 



Hit),piit)Ps 



— I dti Trieads 



to 



(16) 



A significant simplification of (fT6!) is achieved with the so-called Markoff approximation. 
We know that the contacts impose large reservoirs to the SWCNT and assume that they 
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induce fast relaxations inside it. For this reason, the current state of the tube can only 
be influenced by its very latest history, on a timescale not exceeding a certain range r. 
Stating that r is sufficiently small, we may safely average p^{t) over a time period of the 
order of r; this enables us to replace pQiti) in the double commutator in (fT6|) by PQ(t). 
The approximation is also valid when respecting the fact that we examine a static DC 
circuit, where the detailed dynamics on short time intervals need not to be taken into 
account. 

We arrive at the equation of motion we want to work with when we send to —>■ —oo 
(because we are merely interested in the longterm behaviour of the system) and use the 
abbreviation t2 = t — ti : 



h 



Tn 



eads 



HR{t),pQ{t)pSpd 



1 



Tr, 



eads 



dt. 



i^^(t), H^{t-t2),pi{t)psPd 



Removing the commutators, we obtain 

/>Q(t) = Trieads (^AWPoWPsPd + Z^.c.) - (17) 
Trieads^dt2 (^[H^{t)Hi{t - t2)pUt)psPd - HUt)p'Q{t)psPdHUt - ^2)} + h.c}j . 

We can now substitute the explicit form of (given in Eq. ([3D), and of Hj^ (from 
Eq. ([5])) in Eq. ( ITTl) and get rid of the traces. The terms containing Hj^ involve 
no lead operators, such that the trace just makes the product psPd vanish. For the 
contributions with Hj^, the lead operators must be brought next to the density matrix. 
This is done by using the cyclic property of the trace and commuting them past two 
SWCNT operators in each term; this, remarkably, means that it does not play any role 
whether the operators are commuting or anticommuting. 

Then one can introduce the correlation functions 

i^U^^),, ■= Tri-ds [^Ui.,PsPd) , (18) 

and exploit Trieads (^/<7,^Z'<7,,P^Prf) = Trieads {^LM'a,,PsPd^ = 0, 

Trieads (^^laM'a.PsPd) = for lai ^ I'av . (19) 

With the abbreviation £i„^{f,r[t2) 

J^iaXr,r't2) 

Eq. (ITTI) is changed to: 



:= Ur)T;{f') (^^i,^if)^Jf:-t2))^ 
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lai 

la I 

X ( \^£i„Xr,r'M)i!i,X^,t)i!c,aXr:t-t2)p'Q{t) + 

+^,<,,(r,f:t2)^0a,(r,t)^^(f:t-t2)p0(t)} + h.c. - 

+ £l^{r,r'M)i>Q,Xr.t)plS) ^Ui^'.t - t^)] - h.c 



(20) 



Moreover, we are able to write out the time dependences of the operators according to 

independently of all indices but / G {s, d, ©}, because all parts of Hq besides Hj 
commute with 

3.2. Equations of motion in the SWCNT energy basis 

The starting point for this subsection is Eq. (I2U]) . which is used to derive the equations 
of motion for a single element (p0(^))„^ of ths reduced density matrix in the SWCNT 
energy basis. 

At first, we apply two more approximations that will ease our tasks considerably: 

(i) We presume that a SWCNT generally is in a pure charge state, i.e. it is filled with 
a certain number of electrons and thus density matrix elements between states of 
distinct electron numbers are set to zero (we deal with bias voltages not exceeding 
the height of the Coulomb diamonds and permanently measure N in our circuit via 
the gate voltage). 

(ii) Being interested in the static be- 
haviour of our system, we can ne- 
glect fast oscillating terms arising 
from the exponentials in (12T1) for 
/=©. This secular approximation 
completely decouples the time evo- 
lution of matrix elements between 
states degenerate in energy from 
all matrix elements between states 
non-degenerate in energy (fig. [3]). 




Figure 3: The RDM acquires a block 
diagonal form. Light squares: particle 
numbers fixed. Dark squares: also 
degenerate in energy 
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To calculate the current, merely the diagonal elements of the density matrix, the 
populations, are required, and as a consequence it is actually sufficient to deal with 
block matrices Pqe^ni which are restricted to the Hilbert space spanned by states 

] e \EnN) := ||iV') : \N'\ = N, {N'\Hq\N') = En} 

with fixed energy and particle number (dark squares in Fig. [3]). 

Eq. fl20|) then reads, introducing Block- Redfield tensors in order to simplify the notation. 



(PeENNi't)) — Rnm,jj' {p&Enn{^)) jji + ^nm.kk' (^OSfcAf+l (*)) fcfc/ + 

\ / 71771 -^-^ 

jj' kk' 

+ (22) 
a' 

The sums in ( 122|) run over states with fixed particle numbers: J, \E]\fN)}, 
k, k' G {|A^+1)}, i, i' G {|A^-1)}, where it is the secular approximation which additionally 
fixes the energy oi j,j'. 



The Bloch-Redfield tensors are defined as follows: 



^nm,jj' ■ ^ ^ (^^n 



I ie{|7V-i>} 

k G {|Af+l>} 



p(+)Ar7V-l j^(+)7VAr+i 
l{niij) l(nkkj) 



p(-)AriV-l j^{-)NN+l 
l{j'nm) l{j'kkm) 



^nm,ss' ■— 2-^ l{s'mns) ~^ l{s'rnns) J ' y"^"^) 



I 

with the rates 



^ 5^y"dV A^rl sgn(aO (^5,,,(f))^^(^Q.X^l)^,^. • (24) 



We like to point out that up to now, the relations we deduced are very general ones, 
as we did not exploit any property specific for SWCNTs. Eq. (!22|) already shows the 
features that are important for us: the time evolution of the diagonal matrix elements of 
P&EmN will couple to some elements of Pqe^n, PeEiN-i and Pqe^n+i- Which elements 
are actually involved can be mapped out by transforming the expression ( 12^ further. 
First, the explicit form of the correlation functions (fT8|) needs to be determined. This 
and some consecutive steps are carried out in the appendix, and as a result the rates 
are changed to flA.2p : 
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■^{a)NN+l 
l{nkk'j) 



^{a)NN-~l 
l{nii'j) 



X 



r ai 



' nk \ / fc'j 



Di^XEkj) fi{Ekj) + a- 



X 



vr 



7r 



oo e - Efcj 



(25) 



Here, 



Icri I \ 
tot \e) 



e: 



exp 



totle 



+ eVi-E 



Fl 



' (z f) 

is the Fermi function in lead / and has arisen from the correlation functions. Ep j is the 
common Fermi level for the two spin species ai = +/ and ai = — / in contact / without 
any bias voltage applied (for calculations we will assume that source and drain are made 
of the same metal, which allows to shift the energy scale such that E^j^^'^ = sj^^f = 0). 



Further, we have abbreviated Eab = Ea — E^ and $/ , Ri ( lA.ip are real values into which 
the integrations over space, the tunnelling and reflection parameters as well as the wave 
functions from the decompositions (jl]),([9]) of the electron operators have been absorbed. 
To do so, the assumption that both T/(r) and A/(r) are small away from the tunnelling 
contacts was needed. The new parameters and Ri scale with the strength of the 
tunnelling, respectively the reflection. 
Finally, ^ denotes a principal part integration. 

Now let us have a closer look to the rates ( l25i) . As indicated by the notation, r^^"*"^ is 
related to transitions — A^+1; that is why its real part involves the product of the 
density of states Dicri and the Fermi function fi, which refers to the number of electrons 
that can potentially leave the contact. r^^~^ correspondingly contains 1 — which 
accounts for vacancies in the lead. 

The imaginary parts also include a product of Fermi function and density of states, but 
the energy arguments of the functions are integrated over: there is no restriction to 
energetically permitted transitions. That is why we call these virtual; the real parts, 
however, represent processes where the energy is conserved. 



Still, the electron operators in Eqns. ( 1251) are given for spins along the quantisation 
axes of the two distinct lead coordinate systems. To go on, we must reexpress them in 
a common basis. 
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4. Spin quantisation axes 
4.I. N on- collinear magnetisations 

A clever choice for the SWCNT spin quantisation axis Zq is the direction perpendicular 
to the plane spanned by fhs and rhd [21]: using the coordinate systems drawn in Fig. HJ 
we will get particularly nice matrices ^7;,©, which transfer the electron operators to 
their representation "ipl^^ , where the spin is quantised along Zq. 




source contact drain contact 

SWCNT 



Figure 4: The X0-axis is chosen such that it bisects the angle B between fhs and rhd- In each 
case, z is the quantisation axis in the corresponding part of the system. 

To calculate Ui^q, we rewrite the three-dimensional coordinate system basis 
of contact / in the nanotube basis: 

(eOo = Ui^Q {ei)i = Ui^Q (60)0 , 

where e G {x,y,z}, Ui^q G 50(3) and (e)//0 indicates that e is given in the 
lead/SWCNT coordinate representation. 

The second equality makes clear that the coordinate transformation matrix Ui^q is just 
the 3D rotation matrix that rotates the basis {xQ,yQ, Zq} of the nanotube coordinate 
system onto the one of lead I, {xi,yi,zi} . 

Ui^Q itself is simply a product of two rotation matrices of different angles around different 
axes. Carefully considering the directions and the resulting signs for the angles of the 
necessary rotations, we find: 

Ui,e = Us>jei)Ug^{7r/2) with 61: = 

Here Xq denotes the Xq axis, having undergone the first rotation. With some basic 
quantum mechanics we obtain the corresponding transformation matrices for the 
electron operator ipl^^ : 

Ui^Q = (^1 cos(-^;/2) + ia^ sm{~9i/2)) (l cos(45) + iay sin(45)) = 

_J_f +e+i^'/2 +e-'^'/2 
~ 7! [ -e+'^'/2 +e-'^'/2 



-9/2 I 
-e/2 I 



s, 
d. 
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Using ( 1271) it is straightforward to evaluate 
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(27) 



+ 



(28) 



where all uninvolved indices and arguments have been dropped for clarity. Defining 







Z(T0Cr' 



cr«=T,cr^=i and A.^.^ := <( ^^'^^ + ^^'^^ ? ^ !f (29) 



the rates in Eq. fl25|) are reformulated for one last time: 







j^{a)NN+l 
l{nkk'j) 



■^{a)NN-l 
l{nii'j) 



h 



de 



(i?,.) (i-//(i?,.)) 



1 

a — 

TT 



de 



(30) 
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4-2. Collinear magnetisations 

Before we explain how to solve the final master equation, we should have a look at the 
special case of collinear contact configurations, where it is clever to use the common 
quantisation axis of the leads also for the coordinate system inside the SWCNT (fig. [5]). 



(a) 



ys 





(b) 



source contact 



SWCNT 




Figure 5: Considering merely collinear contact magnetisations, it holds advantages just to 
choose one of the lead coordinate systems as SWCNT coordinate system, too. 
(a) depicts the parallel, (b) the antiparallel case. 



The matrices we need for replacing Eq. (12^ when we work with the coordinate systems 
sketched in Fig. [5] are 



Ui,e 



1 cos(O) + ida; sin(O) 
1 cos(— 7r/2) + iax sin(— 7r/2) 



The electron operators transform accordingly. 



1 








1 





— i 


— i 






l = d or 1 = sA9 = 
I = s A 9 = n . 



(31) 



4- 



and instead of (128|) we obtain 



',t 



l = doTl = sAe = 

I = S A 6 = TT , 



(32) 



+ Di_^ljj^.ljji l = doTl = sAe = 

+ l = s A 



TT . 



(33) 



The corresponding rates are not much different from (125|) : 
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■^{a)NN+l 
l{nkk'j) 



r o-Q 



k'j 



Dla^iEkj) fl{Ekj) + 



f — Fy 

CO -'^Jt 



{a)NN-l TtL^ 



l(nii'j) 



h 



ro-Q 



r (T0 



—a — 

7T 



with 7i(6') := Si^sS9,n and 



"4^ fur. 



(34) 



D 



l(7Q 





(70 


= T 




(Tq 


= i 




(T0 


= T 




(70 


= i 



l = d or l = sAe = 
I = s A 9 = n . 



(35) 



The important point is that in Eq. fl33|) . no terms with mixed spins occur at all, which 
will make setting up the necessary equations easier than in the general case for an 
arbitrary 6. In particular, all coherences drop out and for this reason, also the imaginary 
parts the rates fl34p still contain will vanish: in our regime of weak coupling, neither 
the virtual transitions nor the boundary reflections have any influence for coUinear 
configurations. 

Now we can start an analytical examination of the low bias regime, i.e. the linear 
transport. 



5. Linear transport 

What we are going to do analytically is to consider linear transport, where eVJ <^ ksT , 
at thermal energies ksT -C eo < ejn ^jn. In this case, both fermionic ([Jj]) and bosonic 
excitations can only be created virtually. That is why \N,m) = \N, 0) is certainly true 
for all states contributing to transport in our limit, and caring only about different 
fermionic configurations we symbolise the associated eigenstates by |A^). Due to the 
absence of any excitations, the four bands of the SWCNT have to be populated as 
equally as possible, i.e. their electron fillings , A''^^ , A''^^ , A''^^ may at most differ 

II In the presence of fermionic excitations the four tctq bands are no longer filled as equally as possible: 
there are at least two band whose electron fillings differ by more than one. Still, however, it is the 
lowermost states of each band which must be populated. "Holes" can solely be produced by bosonic 
excitations 
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by one. 

Involved states |A^) are thus for fixed |A^| characterised by the energy 



E 



:= min{E^ : |A^| = A^} . 



N 



They are the groundstates, i.e. the states with lowest possible energy for a fixed tube 
filling, and therefore we call them lowest energy (ground) states (LEGs). For a SWCNT 
without a band offset 5, Fig. [6] shows all these LEGs which can, under the given 
conditions, a priori be involved in transport. We will soon explain that also for the 
virtual processes, only transitions between the LEGS contribute in the end. 
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Figure 6: All possible lowest energy (ground) states for a SWCNT without band offset, filled 
with 4N + A electrons. 

At low bias voltages, the current flow across the SWCNT is Coulomb blocked, unless the 
applied gate voltage Vgate aligns the states \E^p N) and \E^p^-^^N+l): resonant tunnelling 
then permits to add or remove the 1st inside the SWCNT without any energy cost. 

In this article we will just consider the situation 5 = as shown in Fig. [6l because this 
is sufficient for describing a SWCNT also in the case of a flnite band offset 6 > ksT 
completely. Namely, if the thermal energy cannot overcome the energy gap between 
the R- and the L-band, resonant tunnelling can solely be allowed for one of those two 
bands: at each value of the gate voltage, we just deal with a speciflc of two independent 
single level quantum dots, while the other one is not of importance. Any regime of such 
a system can be mapped on either the resonance A = 0<->-A = lorA = 3^A = 
for a zero band mismatch. 

To calculate the current, we need the populations of the LEGs, 

^6{|<>7V>} 
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which represent the probabihties to find the SWCNT filled with electrons. As we want 
to examine the DC long term behaviour of the system, we claim that the populations 
do not change with time. If the value of Vgate is tuned to admit resonant tunnelling, this 
stationarity demands: 

l=s,d l=s,d /otN 

= = - 5^ + 5^ ^r^+^p^ • 

l=s,d l=s,d 

Here, j ^A^+i->Af (jg^ote the rates by which electrons tunnel into (A^ — *• 

A^ + 1) / out of (A^ + 1 — > A^) the nanotube at contact /. 

The question is now to which elements of pQ^it), pQ^^i(t), Pq7v_i(^) the time evolutions 
of the populations, Pjy and Pn+i, couple. For an answer, the rates - either Eq. fl30|) or 
( !34|) - must be consulted: we find that their real parts always include Fermi functions 
with an energy difference as an argument. While Vgate is chosen such that 

all other energy differences are at least of the order of the charging energy Ec and due to 
the condition ksT <^ eo ~ E^., the real parts of the corresponding rates are suppressed 
enormously. Not the rates, but the Bloch-Redfield coefficients fl23l) enter the equation 
of motion for the density matrix. By its definition, i^^^^i is real and thus the only 
nonvanishing contribution is the one of R^mtk' with k, k' being LEGs: merely elements 
of the density matrix P^^(o) ^_,_]^(^) involved. For pQ^{t), the Kronecker deltas in 

Rnmjj' demand either n = j or m = j'. Hence, En = Em = E^^^ and consequently 
Ej = Ej = E^^^: no other element than the ones of /3^^(o)^(t) itself are participating. 
With the help of Fig. [6], finally the concerned Bloch-Redfield coefficients can be figured 
out. Before, however, the role of the electron operators must be revealed. From [8], it 
can be extracted that: 



If nonzero, the sign of (1381) depends on whether an even (+) or an odd (-) number of 
electrons with the same energy sits in the bands prior to ra©. To define prior, one has 
to choose an arbitrary but fixed ordering of the bands, e.g. L be prior to R and r | be 
prior to f f. 

For our case, a, c G {\E^^'^ N)} in Eq. ([38]). What about h E {|A^+1)}? Still there is a 
sum over all TnUm and Tnkkm in Rn^jj' ' with i e {\N-1) and Tnkkm with A; G |A^+1). 
For the real parts of the rates, the energy argument of the Fermi functions restricts 
these i and k to LEGs, but for the imaginary parts, this energy argument is integrated 
over and so there is no forbiddance for fermionic or bosonic excitations a priori. After 
some closer examination, however, it turns out that really all non-LEG contributions 
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of the four rates Rnmjj' contains drop out exactly. Especially, one has to comprehend 
that if for example in 

^ J jk \ ^'^0/ kj' 

k G {A^+1)} was a non-LEG while j and j' are LEGs, the expression can be nonzero 
only for raQ = r'a^: if we create an electron to reach an excited state and wish to 
return to a LEG, there is actually no other possibility than to remove again the electron 
we added. Of course, j = j' then, such that the concerned rates ^[^^kkj)'^^ ^\{fkkm}^ 
in Eq. fl23l) become '^[^nkkn)^^ ^^'^ '^[imktm)^ '■• imaginary parts of the latter, however, 
cancel each other. 



After all, we can set up Eqns. ( 1371) for the four different tunnelling regimes 
A ^ A + 1, where A := A^mod4. Once we have done so, it is easy to extract the 
current by just omitting the sum over one of the leads: 

|/aa+i| = e |Af-^+iPA - Af+i-^PA+il = e {A^-'^^'Pa - A^^'-^P^^,\ . (39) 

Notice that from here on, we replace the indices by A, because for all the properties 
we will study just mod 4 = A and not the actual value of A^ matters. 

For the calculations we state 

+ Di-iifJ'A) =■ A(/iA) = Di ; 

this is justified since the densities have to be taken at the Fermi edge Ep^i of the lead 
metal and Epj ^ /iA, so that we can assume them to be approximately constant within 
the range we consider. 

Moreover we define the contact polarisation: 
5.1. Resonant tunnelling regime A = ^ A = 1 

If A = 0, there is only one LEG, with all bands equally occupied (see Fig. [6]). We have 

pLoW = (pLoW) 11=^^0, 

which is the probability to find the SWCNT filled with A^ = 4iV , N e N electrons. 
Note that for the new variables we drop the argument (t) to save some space in the 
following. 

For A = 1, things are a bit more complicated, because the LEG is fourfold degenerate: 
the excess electron can be placed in each of the four bands L ^ , L I, R ^ , R I ; 
consequently, the density matrix consists of sixteen elements. Eight of them, however, 
can immediately be set to zero, as we know that our system is unpolarised with respect 



Spin transport across carbon nanotube quantum dots 



21 




to L- and i?-band. This forbids any transitions and thus coherences between those two 
bands: 

pi=iW = 

Here, due to the indistinguishabihty of L- and i?-band, we are allowed to set 

[ P^^^ p(i)e*"''' 

^A=1,lW = ^A=l,ij(^) = I p{l)Q-iaW p{l) 

where P^^^ , Pj^^ are the probabilities to find the SWCNT in a single-electron spin-down 
(I) respectively spin-up (t) state. Correspondingly, 

Pi := 2P|^^ + 2P^^^^ 

is the total occupation probability for one electron. 

Furthermore, the density matrix is hermitian and that is why we could define 



for the off-diagonal elements. 

The meaning of these quantities is revealed when we extract the information about the 
average spin S*^^^ on the quantum dot from p^^^; a single spin-1/2 particle that can 
either be in the spin-up or the spin-down state would be described by a 2 x 2 density 
matrix and the average spin would be given by a trace with the Pauli matrices. Now 
we have additionally the two f-bands, but those decouple completely, such that we can 
obtain the components of the average spin by 

where aj are the Pauli matrices, and therefore 

^(1) = 1 (2j9(V"^% 2pWe--'^') = 2p' 

^(1) = 1 (2p(i)e'"'^' - 2p(^)e-'"^^') = -2p 

1 
2 

If we chose our spin quantisation axis inside the SWCNT such that ZqWS'-^K then 
s!i'^ and Sy^^ and with them all off-diagonal elements of pa=i would vanish: we had a 
coordinate system where the hermitian density matrix is diagonalised. For arbitrary 6 , 
we are not able to give the direction of a priori and thus cannot make any use of this 
insight. For ^ = and 6 = tt, however, we have already found the right axis and could 
obtain diagonal equations when using the rates instead of (130|) . But actually the 
A = ^ A = 1 regime is analytically accessible for any 6 and as it is instructive, we 



W cos(a«). 



sin(a(i)) , 



5(1) = - (2P(1) - 2PW ) = Pf ) - Pf ) 
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want to set up the general equations, in terms of the physical quantities Pq, Pi, S^, Sy 
and Sz, which are five independent variables determining the density matrices Pa=o and 
Pa=i • 

The equations for Pq and Pi are not independent (Pq = —Pi); but together with the 
normalisation condition 

Po + Pi = l (40) 

we have the necessary set of five equations, where the physical meaning of the single 
terms becomes obvious: 



^1 = ¥ E ^^^^ Ufiif^A=i)Po - (1 - /K/iA=i)) Pi - 2Vi (1 - fiifiA=i)) (5^(1) 



d 



nil 



electrons electrons 



,,. . , ... , difference m chemical potential 

tunnelling in tunnelling out ^ . , ■ , 

tor spin— up and spin— down 



dt h 



l=s,d '"v,.^^^ 



2fi{fiA=i)Po - ^ (1 - /Km=i)) Pi ) rhi 



spin accumulation 



- (1 - //(/iA=l)) 5^(1) - ^ %{fiA=l, /iA=2) C 
^ ■' . TT V 



nil X 



5?(i) 



spin relaxation ^ 

spin precession 



(41) 



The probability Pi of the SWCNT to be occupied with AN + 1 electrons grows with the 
rate of electrons tunnelling into the tube already containing 4N electrons and decays 
with electrons leaving the quantum dot occupied with 4A^ + 1. Additionally, the average 
spin of the electrons inside the SWCNT interacts with the magnetic fields in the contacts, 
which yields a special term accounting for the difference in the chemical potential for 
spin-up and spin-down electrons. 

That is why we have also to solve the second equation for S^^\ The time evolution of the 
average spin is affected by three contributions: the net number of particles accumulating 
on the tube brings a spin polarised along the lead magnetisations with it; electrons 
tunnelling out take some spin with them such that the spin inside the tube relaxes, and 
finally the electrons feel the magnetic field of the contacts, what makes them precess a 
httle bit. 

Besides of a factor 2 (due to the fact that we have the distinct r bands) and the precise 
form of ^i, these are the very same equations obtained in [21] for a single level quantum 
dot. The function 

J Ve-/^A=i e-AtA=2/ ViDi 
summarises all the imaginary contributions the rates contain; this means that it merges 
our two exchange effects: the interface backscattering processes - which have not been 
included in [2T] - and the virtual transitions. The exchange is responsible for the 
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precession of the total spin on the quantum dot. 

The principal part integration appearing in can be evaluated by a trick: 
for a Lorentzian-like shaped function L(e) = L — — 



2_i_e;2 , 



E max {ksT, \ E\) 

For By/ oo, L(e) is constant and approximately we may set our function Di^\^{e) = 
Di+^{e) — Dl_^{e) = ViDi constant under the principal part integral. This yields for 

cnf ^ 1 max (/csT, I /iA=2 1) , Ri 

max(fcBT, |/iA=i|) ^iViDi 
' V ' 

Why do we employ this estimation for the digamma function [20]? In our analytical 
results, we will find accompanied by a Fermi function, which at low temperatures 
dominates the gate voltage evolution of the product. Therefore, we will, though (H2l) 
cuts the peaks of the digamma function, without any exception obtain smooth curves 
(which indicates that we make no noticeable error). 

For the simplest case of identical source and drain tunnelling contacts, 

D, = Dd=:D,Vs = Vd=:V, = =: $ , R, = =: R 

-and thus = =: we can give a nice analytical expression for the low 
bias current in dependence on both the magnetisation angle 6 and the bias voltage 
Vgate- The zero point of the latter we fix for each resonance A ^ A + 1 such that 

V^ate = -v^ -^A+l ~ -^a'' ~ ' "which means Vgate = —j^A+i- 

/oi(^, V..... Hias) = /oi 1 - ^2(,..,,;i; , Hias , (43a) 

where Vd = —Vg = and 

hkBT 1 + 3/(/iA=i) 

just as it was derived in [8] for a SWCNT attached to unpolarised leads. Fig. [7^ shows 
Joi in dependence on V^ate around its resonance for several 6 and the corresponding 
TMR curves 

TMR (9) = ^Sate) - /aA+i(6', Vgate) ^^^^ 

-^AA+l(0, V'gate) 



at A = 0. Notice that here we have defined a TMR for arbitrary angles. Eq. ( l44l) 
expands the definition we gave in the introduction for the special case of 6' = vr. 
In Fig. [Hk, V'gate is fixed at three different values and the normalised current /oi(6')//oi(0) 
is plotted against 9. The parameters we employ are given beneath the plots. 
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Figs. U}p and [Hb belong to the A = 3 ^ A = resonance and are mirror-symmetric to 
[7^ andlHK with respect to V^ate = 0. This fact results from the symmetry in the LEGs: 
for A = 1, there is one excess electron that can be put in any of the four bands, and for 
A = 3 it is just the same situation with one hole instead of an electron. The resulting 
formula for the current is actually: 



-^30 ^gate, Hia 



/30 I 1 



P^sin^(l) 



(I) 



Vt 



bias ' 



and 



kksT 1 + 3/(-/iA=o) 



(45 a) 



(45 &) 



Unless specified differently, we will employ for these and all following plots the 
parameters listed below in Tab. [TJ 



Polarisation 


V = 


0.6 


Reflection 


R/{^VD) = 0.1<Po 


Tube length 


Lt = 


580 nm 


Temperature 


T = 


20 mK 


Thermal energy 


ksT 


^ 1.73 fieV 


SWCNT charging energy 


Ec = 


-- 9.49 meV 


SWCNT level spacing 


eo = 


2.89 meV 



Table 1: Parameters used for the plots (if not specified differently) 




Figure 7: Currents /oi(^ate) and /3o(V"gate) for different values of 6, together with the TMR 
(j44p . The symbols represent numerical data, which is perfectly fit by the analytical lines. The 
corresponding angle-dependent TMR curves are plotted as lines as well. 
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Figure 8: Normalised currents /oi(6')//oi(0) and hoiO) / l3o{0) for different values of Fgate- 



Besides for the mirror-symmetry with respect to V^ate = 0, there is only one more 
difference between (143 ap and (145 ap : the argument of ^ is /iA=2 = — AiA=i for the 
A = ^ A = 1 regime, but atA = 3^A = 0, additionally the level spacing comes 
in: /iA=i = Ec + eQ — /iA=i- The resulting deviations, however, are washed out by the 
squares of the trigonometric functions and the polarisation in our formula and thus are 
not noticeable in figs. [7] and [HI As already mentioned, the dominating gate voltage 
dependence of the Fermi function by which ^ is divided regulates many features of the 
plotted curves. For instance, it determines the roaming of the peak maxima towards 
V^ate = for 6 ^ {0, tt}. In the case of collinear magnetisations {6 G {0, vr}) , the maxima 
positions, which are identical, can be calculated easily by differentiating fl43 b\i . They lie 
a bit off-resonance, at 

_ f-A;Brin2 ^ -1.2/ieV for A = ^ A = 1 , 
ea gate - j ^^^^1^2 ^ +1.2 /leV for A = 3 ^ A = . 

Plugging the currents into Eq. (jBl), the evolution of the TMR can be explained. For 
^ = vr , it is constant at the value 1 — V"^. For non-collinear contact magnetisations, 
again the Fermi function in question is the decisive factor: the larger the inverse of 
its square becomes, the closer gets Iaa+i{0) to /aa+i(0) and the TMR vanishes. In 
a physical sense, we can imagine the following picture: below the A = A = 1 
resonance, a spin-polarised electron which is transferred across the quantum dot will 
not stay inside the SWCNT too long and thus not have time to equilibrate its spin, 
which hinders the transport, yielding a nonzero TMR. Above the resonance, the tube is 
mostly populated with 4N + 1 electrons, such that for all but exactly antiparallel config- 
urations, the spin of the excess electron will have equilibrated before it tunnels out, and 
consequently the non-collinear TMR decreases around the resonance with growing gate 
voltage. In analogue, the mirror-symmetric behaviour of the A = 3 ^ A = regime 
can be understood, when thinking in terms of holes rather than electrons. 
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This also explains the dependence of the normalised currents on 6, fig El Transport is 
blocked by an antiparallel polarisation of the contacts, and that is why the plots exhibit 
a dip at 6^ = TT. The more likely a spin-equilibration, the closer is /aa+i(^) to /aa+i(0), 
i.e. the narrower get the curves. So for A = 0, the width of the dip must shrink when 
raising the gate voltage, while for A = 3 it is just the other way round. 



E 




^ TT 

Magnetisation angle 6 [rad] 

eaVgate = 0/ieV R/{'S>VD) = 
— 
iPo — 
4q3o — 

loqio — 



Figure 9: In the case of symmetric 
tunnelling contacts, the dip width of 
the normalised current Jqi (0) / Iqi (0) 
shrinks with growing R. 



We see from Eq. (1421) that ^ grows 
(and hence the dips widths would shrink) 
with R: the refiection processes fortify 
the equilibration of spins. As a result, 
the width of the normalised current curves 
Jqi (6*) /Jqi (0) must shrink with increasing R, 
which can clearly be seen in Fig. [9l 
One should mention at this point that 
we chose the value of R (Tab. [1]) small 
enough to see the differences for distinct 
angles in figs. [7] and [HI For real systems, 
the actual value of R might very well 
depend on how the different domains of 
the ferromagnetic contacts couple to the 
graphene sublattices and probably R differs 
from sample to sample. At the current state 
of the experimental art, where adjusting 
a well-defined angle between the lead 
magnetisations is still a big challenge, we 
find it important to visualise the qualitative 
impact of the different quantities. 



We can furthermore give the analytical 
expressions for the probabilities Pq 
and Pi to find the SWCNT populated 
with N = AN respectively N = AN + 1 
electrons (fig. fTOl) : 



Pi 



/(-/iA=l) 
l + 3/(/iA=l^ 

4/(/XA=l) 
l + 3/(/iA=l^ 



(46) 




N = AN 

N = 4:N+1 



-12 -6 6 12 
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Figure 10: Populations Pq and Pi around 
the resonance A = <-> A = 1 



Of course, the occupation probability for an additional electron grows with the gate volt- 
age and it is worth to be stressed that the Eqns. ( l46l) do not depend on anything else 
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than Vgate- The populations become Pq = -Pi = | at eaV^ate = —ksT^, where also the 
maximum of the resonance is located. 

The components of the average spin in the A = 1 state follow simple dependences as 
well: 



-1 



/q(1)\ e p -D • ^ , ^2 'P^(/^A=1,M=2)C0S^ ^ 

{Sy) = -1— qKias PiVsm-\l+V ^ ^ — 

_^ y(/^A=l,/iA=2)cOsf 

In Fig. [niwe show (4^^) and i^Sf^) for legate = and legate = ±5.2 /ieV. 



eaVgate = — 5.2/ieV eaV^ate = /ieV eaV^ate = +5.2 /ieV 




< 0f7r^27r 0f7r^27r 0f7r^27r 



Magnetisation angle 61 [rad] 

Figure 11: The peaks of (5^^^) and {si^^) get more and more pronounced when raising the 
gate voltage. 

The total spin grows with the polarisation V and linearly with both the bias voltage 
Vbias and the occupation probability Pi. Due to the choice for our coordinate system 
(remember Fig. |4]), the X0-component of the average spin is zero. The fact that (4^^) 
peaks at 6' = vr is clear, because the i/Q-components of the lead magnetisations are op- 
posite and scaling with sin ^. It is the spin precession, which tilts the accumulated spin 
out of plane and therewith gives rise to its nonzero ^-component; consequently (Si^^) 
is proportional to its sign changes with the one of the cross product in (j^Tl) . at 6* = tt . 
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What happens if the transparencies of the two contacts differ? Let us set 

Ds = Dd = D, Vs = Vd = V, = = . 

Only for the parallel and the antiparallel current, a nice analytical expression can be 
given: 

loi{0 = 0, Vgate, e) = Ioi{9 = 0, Vgate) , 

T (ft = ^V .\= 2eil + e)Ue = n,V,,,,) ^^^^ 

' ' ((1 -V)+ e{l + V)){{l + V)+e{l- V)) ' 

Obviously, the ratio j°^|g~Q'g| is again constant with V^ate and the positions of the peak 
maxima are still the same as in the symmetric case. 

Fig. [T2b presents the numerical results for the normalised current at two gate voltages 
(±5.2yU,eV), for two strongly differing values of R (for simplicity we assumed that the 
reflection scales with the transparency, i.e. Rs/^s = Rd/^d =■ R/^)- 
We see that the four curves can hardly be distinguished from each other, and especially 
Ioi{0°,e) > loi{0,e) MO G [0°,360°] : there cannot appear any negative TMR. Fig. [T2b 
shows the analytical peaks for ^ = 0° and 9 = 180°, along with the numerical data, 
which additionally gives the peak for 6 = 140°. /oi(0°,£:) and Joi(180°,e) will, as Eq. 
(148|) tells us, not differ for distinct R. Also for the non-coUinear magnetisation, altering 
R over four orders of magnitude obviously has no notable effect. No negative TMR can 
be found throughout the full range of the gate voltage, and in particular, the dependence 
of the non-coUinear TMR on the gate voltage is even suppressed in comparison with the 
symmetric coupling. 
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Figure 12: (a) /oi(6')//oi(0) and (b) /oi(V"gate) together with the related TMR for = 0.3 ^-rf. 
The normalised current is plotted at different gate voltages V^ate and hugely varying reflection 
parameters R to show that neither V^ate nor R influence its angular evolution much. 
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Figure 13: (a) /i2(^)/A2(0) and (b) /i2(V^gate) together with the related TMR around the 
resonance A = 1 <-> A = 2. Notice that in contrast to the A = A = 1 transition, the 
TMR exhibits a maximum in the vicinity of the conductance maximum. 

Because of the large degeneracy of the involved states, the A = l^A = 2is not even for 
the symmetric case accessible analytically if 9 is arbitrary. We therefore show numerical 
data for 9 ^ {0, tt}. The regime A = 2 ^ A = 3 has again a mirror-symmetry to the 
one we treat here and therefore omitted. For the P and the AP current, an analytical 
solution can be provided: 

T (nr^v ^ _ 37r'e^ /(/iA=2)/(-/iA=2) 

hkBT 2 + f{nA=2) (49) 

h2{9 = TT, Vgate) = (l " V^) /i2(^ = 0, legate) , 

in the symmetrical case. 

The TMR at 6' = vr is consequently constant, at the same value as in the A = ^ A = 
1/A = 3«-^A = regimes. The curves for 6* 7^ vr, however, qualitatively differ, due 
to degeneracy in the f bands, which allows a nonzero average spin in the A = 2 state. 
Below the resonance, there is mainly one single excess electron inside the SWCNT, 
with equilibrated spin. During charge transfer, an additional, spin-polarised electron is 
entering, but it can be the first electron (the one with relaxed spin) which leaves the 
tube. That is why polarised non-coUinear transport is easier, and hence the TMR much 
smaller, than below the A = ^ A = 1 peak. Approaching the resonance, however, 
the time intervals between the described process and the next charge transfer shortens 
so that the A = 1 state does not persist long enough to allow a sufficient equilibration 
of the excess electron's spin. Therefore the TMR increases. Beyond the resonance, the 
A = 2 is stable long enough let the total spin of the two electrons equilibrate. Both 
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excess electrons can now be involved in polarised transport and so the TMR drops to a 
value lower than below the resonance. Still it stays nonzero, as an outtunnelling particle 
always leaves a single unequilibrated spin inside the SWCNT. 

Again, the current peaks are located off-resonance, at 

In 1.5 

eaVgate = —ksT — - — ~ —0.35 fieV (50) 

for 6 G {0,7r}, as the derivative of fl49p reveals. For angles in between, the peak 
maximum again moves slightly towards eaVgate = 0. 

The normalised current luiO) in Fig. [T3](a) looks rather familiar at eal^ate = 0, but 
nevertheless does not follow the analytical law we found for Ioi{9). Off-resonance, Ii2{0) 
obviously no longer exhibits any handsome dependence and as explained, the polarised 
non-collinear transport is hindered at the resonance so that the width of the dip is for 
eaVgate = fieV larger than for both eaV^ate = ±5.2 /ieV. 

For an asymmetric setup, one finds for varying R as less deviations as in the A = ^ 
A = 1 regime. Again, the current in the nearly antiparallel setup with 6 = 170° hardly 
differs from the one for an exact AP configuration with 6 = 180° and the gate voltage 
dependence of the TMR is suppressed. Still, /i2(0°,£:) > Ji2(170°,e) is guaranteed; 
the equations for the P and AP configurations can again be solved analytically and we 
recover our previous finding: 

Wi(g = ^) {l-V'){l+e)' 

/aa+i(^ = 0) {{i-V)+e{l + V)){{l + V)+e{l-V))- ^ ' 

So the ratio of P to AP current is independent of the gate voltage and constant even 

for a possible band offset. 



6. Nonlinear transport 

Being interested in the nonlinear transport behaviour of the ferromagnetically contacted 
SWCNT quantum dot, one has to take into account the various non-LEG states of the 
system, together with all the arising coherences. This makes an analytical calculation 
impossible, but the numerics delivers the relevant data respecting the influence of both 
real and virtual transitions with the energy 3eo of up to three (bosonic or fermionic) 
neutral excitations. 

An inclusion of higher excited states would multiply the computational cost, but 
anyway, their contribution for bias voltages below the second excitation is rather small: 
ksT <^ eo, so that thermal excitation is inhibited and the probability of bias voltage 
driven multiple subsequent excitations naturally decreases with the amount of required 
subsequent excitations [8]. 

Again, we can confine ourselves to distinguishing between the two regimes A = 
A = 1 and A = 1 A = 2, because the symmetry arguments still hold for all states. 
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Figure 14: Nonlinear current Iqi{9 = 0) below the first excitation. For the computation of 
this plot, transitions with a maximum energy eo were taken into account. 



Fig. [13] shows the current for bias voltages eVi <^ eq , i.e. far below the first excitation. 
We can nicely see how the Coulomb blockade diamonds for A = and A = 1 emerge; 
the higher the bias, the wider gets the range of the gate voltage within which transport 
is allowed. This brings about the typical rhomb shapes. 
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Figure 15: TMR vs. gate and bias voltage for = 140° and 9 = 180° around the 
A = A = 1 resonance. As we have found it in the linear bias regime, Fig. [71 the non- 
collinear TMR {9 = 140°) drops from a constant value to zero after the resonance, following 
the current. For 9 = 180°, the TMR is practically symmetric with respect to the resonant 
tunnelling regime, where it acquires an increased value. 



Figs. [T^ and [TH ) are colourmap TMR plots for the polarisation angles 9 = 140° and 
6 = 180° in the same regime Fig. [14] covers. In Fig. [T5h we regain what we already 
knew from Fig. [7] for the non-collinear magnetisations, namely that the TMR drops 
from a constant value to zero after the resonant tunnelling regime. Fig. [T5b . however, 
shows that in the antiparallel configuration, for higher bias voltages the TMR is not 
everywhere constant with the gate voltage but increases outside the Coulomb diamonds. 
Nevertheless, just as for the low bias, TMR(6' = 140°) <TMR(6' = 180°). Moreover, 
both cases lack any negative TMR and actually the TMR gets widely constant for high 
bias voltages, as the nonzero currents in between the Coulomb diamonds acquire a 
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constant value. That is why for eVi > e^, we only need to plot the current at one special 
gate voltage, e.g. Vgate = (fig. [16]). 
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Figure 16: Nonlinear current Ioi{9) for 6 G {0°, 140°, 180°} and bias voltages exceeding the 
first neutral excitation. The numerics respected a maximum transition energy of Seo- For a 
non-collinear lead magnetisation, a negative differential conductance appears below the first 
excitation. 

Within the range of the bias voltage fig. [16] shows, the curves for 6 G {0°, 180°} reach 
three different plateaus: below the first excitation, there is a constant current outside 
the Coulomb blockade diamonds. As soon as the bias voltage is high enough to excite 
a neutral mode, more states can contribute to transport and the current jumps to a 
higher value. It is important that from here on, successively by picking up the energy 
the bias voltage provides, multiple excited states can be generated. The probability 
to fall back to a less excited state, however, is always larger than the one for creating 
another additional excitation and therefore it is certainly valid to at an energy of 3eo. 
Indeed, another slight enlargement takes place at that value of the bias, which provides 
enough energy to add an additional electron to a state |A^+1, 1); this is another allowed 
transition above the first neutral excitation. 

The current at non-collinear magnetisations exhibits a further feature: we find a negative 
differential conductance (NDC) below the first bosonic excitation. This behaviour 
becomes more evident for high polarisation, but can still be seen for our passably 
realistic case oi V = 0.6. The explanation for the occurrence of the NDC is the 
decaying infiuence of the virtual transitions with growing bias voltage (due to the energy 
arguments appearing under the principal part integral, see Eq. ( l25i) ). We learned that 
the principal part terms narrow the Im{9) curve (fig. [8]), which means that the heavier 
their infiuence, the close comes the non-collinear current for a certain magnetisation 
angle 9 to the maximum current /oi(0°). Thus it is clear that Joi(140°) approaches the 
minimal current /oi(180°) for higher bias voltages where the principal part is more and 
more suppressed. Notice that here R has an infiuence, but as it does not depend on any 
external voltage, it tends to wipe out the effect. 
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Figure 17: Nonlinear current Ii2{0 = 0) below the first excitation, transitions with a maximum 
energy eg taken into account. 

At the A = 1 ^ A = 2 resonance, most statements from the previous subsection hold 
true as well. Due to the higher degeneracy of the states involved in transport, the 
plateau the current reaches is 50% higher than for the A = ^ A = 1 resonance [8]. 

^ (a) TMR(6i = 140°) (b) TMR (6i = 180°) 

T 

^ -12 -6 6 12 -12 -6 6 12 

Gate voltage eV^ate [/^eV] Gate voltage eV^ate [m^V] 

Figure 18: TMR vs. gate and bias voltage around the A = 1 ^ A = 2 resonance. Again, the 
non-collinear TMR for the polarisation angles 9 = 140° is similar to the one in the linear bias 
regime: before dropping to zero after the resonance, it shows a maximum near the conductance 
maximum. Just as the latter one, the TMR maximum widens with the bias voltage. The 
collinear TMR is qualitatively equivalent to the A = A = 1 case of Fig. [151 but the 
increment of the TMR is for A = 1 A = 2 considerably less, namely about 10%, compared 
to some 40% for Fig. (TJ] 




For the colourmap TMR plots of Fig. [18], the same scale as for Fig. [15] was applied. 
Again, at ^ = 140° we find an evolution of the TMR we expected from the linear 
regime. Fig. [T31 and for 9 = 180° we obtain a qualitatively similar picture as Fig. \T5b 
for A = ^ A = 1. Figs. [T5b and [T8b. however, differ in quantity: the maximum 
value (~ 0.39) of the TMR atA = l<^A = 2is close to the value 0.36 for the linear 
bias, while atA = O^A = lit still rises with the bias voltage from 0.36 to about 0.5. 
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For the higher bias vohage regime, Fig. [T^l we can observe the same facts described 
in section 16. 1[ It is worth mentioning that beyond the first bosonic excitation, the 
differences in the heights of the currents at the two distinct resonances get much smaller, 
because the number of involved states multiplies in both regimes. 
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Figure 19: Nonlinear current luiO) for 9 £ {0°, 140°, 180°} and bias voltages exceeding the 
first bosonic excitation. The numerics respected a maximum transition energy of Seo- 



7. Main results 



Within the framework of our master equation approach, we could deduce equations 
for spin-dependent transport across carbon nanotube quantum dots. In the tunnelling 
regimes belonging to tube fillings AN ^ 4A^+1 and 4A^+3 ^ 4(A^+1), the system 
behaves equivalent to a single level quantum dot. The resonances 4iV+l <-> 4iV+2 and 
4N+2 ^ 4A^+3 show similar characteristics; however, due to the larger degeneracy, 
the currents are increased and the fact that the dot is always populated by at least 
one electron brings about a more complex evolution of the TMR. In our perturbation 
approach (low transmission, which involes also only a small spin-splitting produced by 
boundary reflections), neither for symmetric nor for antisymmetric coupling to the leads, 
a negative TMR can occur. 
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8. Conclusion 

We have investigated spin-dependent transport in SWCNT quantum dots. A spin- 
dependent equation for the reduced density matrix of a SWCNT weakly contacted to 
ferromagnetic leads of arbitrary magnetisations was presented. We demonstrated that 
the SWCNT behaves as a spin-valve single electron transistor and showed analytical 
and numerical results for the current flow. 

Because of the fourfold periodicity for the electron number A^, it is sufficient to 
discriminate between tube fillings with different values of A = mod 4 and, due to 
mirror-symmetries in the SWCNT eigenstates, we could even restrict our examinations 
to the tunnelling regimes A = ^ A = 1 and A = 1 ^ A = 2. 

The analytical analysis in the case of symmetric coupling to the leads resulted for the 
A = A = 1 resonance in an equivalent formula for the angular dependence of 
/oi(6') as [21] and [23] obtained, for a single level quantum dot and a metallic island 
respectively. The total current, due to the existence of the degenerate left and a right 
mover bands, is twice as large as for a single level quantum dot. The maxima of the 
current peaks lie slightly off-resonance; the positions for a parallel and an antiparallel 
magnetisation are identical, but on the way from 6' = 0° to ^ = 180° the maximum 
moves a bit, following a curve bent towards the resonance gate voltage (in the /oi(^ate) 
diagram). The TMR for 6 ^ 180° changes around the resonance: it smoothly drops 
from a constant value, for gate voltages below the resonance, to zero. 
We additionally gave the average spin on our quantum dot SWCNT in the A = 1 state 
and the two occupation probabilities Pa=o and Pa=i, where we find that the latter solely 
depend on Vgate- For a nonzero band offset < ksT ^ 5 we can apply these results, 
because it makes the SWCNT at all gate voltages equivalent to a single level quantum 
dot. 

For 5 <^ ksTi^ eo), the resonance regimes A = l^A = 2is more complex. The 
TMR around the resonance not simply decays monotonously from one constant value 
to another, but shows a peak before decreasing. Nevertheless, all TMR curves are 
strictly positive and also for an asymmetric coupling to the leads, numerical results in 
the different regimes at ^ = 140° reveal that a non-collinear contact polarisation alone 
cannot produce a negative TMR. 

Specifically, we could deduce the general law (jSTj) for the TMR in the case of collinear 
(P-AP) magnetisations. It shows that under strict lowest order perturbation treatment 
of both the tunnelling and the refiection parameters, the linear bias P-AP current is 
even for an arbitrarily asymmetric coupling independent of the gate voltage. In order to 
reproduce a negative TMR as observed by [IT], a spin- dependent energy shift (which can 
be obtained from a non-perturbative treatment of H^) is necessary. Actually, exchange 
effects emerging from a distinguishability of inter- and intra-lattice interactions can be 
source of an intrinsic spin-dependent energy shift, as measurements shown in [16] on 
unpolarised small-diameter tubes and recent theoretical investigations exhibit [29]. The 
feature, however, is only present at the resonances involving A = 2 and therefore not 
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responsible for the required gate-voltage independent spin-splitting. 

For the nonlinear bias voltage regime of the quantum dot SWCNT, the numerical data 
revealed that the TMR for P-AP configurations is no longer strictly constant, but rises 
inside the resonant tunnelling regimes, whereas atA = O^A = l, the effect is much 
more pronounced than at A = 1 ^ A = 2. The non-collinear TMR at ^ = 140° is 
similar to the linear bias TMR, but now it becomes obvious that the changes in value 
take place at the edges of the Coulomb diamonds. 

Tracing the current at the resonances to bias voltages exceeding the energy of the 
first possible excitation, we find for both tunnelling regimes a qualitatively equivalent 
dependence. Besides for the large jump at the first excitation, another small one can 
be found when the bias voltage reaches a value that provides the energy to enable 
transitions to states with a second additional electron. A special feature of the non- 
collinear polarisation is that a negative differential conductance appears below the first 
resonance due to the decaying infiuence of virtual processes. 

* * * 
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Decomposition of the correlation functions 



To rewrite the rates ( l24l) . 



nk\ / fc'j Jo 



we first have to determine the explicit form of the correlation functions (fT8|) : 
and in analogue 



fi{E[:i le) = 1 + exp ^^^1^ ^ 



Here, 



is the Fermi function in lead /, where E'^pP is the common Fermi level for the two spin 
species ai = +/ and ai = — ; in contact / without any bias voltage applied. 
Additionally inserting the decomposition of the SWCNT electron operator, Eq. (Q, and 
introducing the quantities 

^irAe):=Y,^gniFF') /"dV/"dV V{r^%ir')Y,'^>h(^Mn^HniF)r]Fif^^^^^^^^^ 

FF' J J 

and Rirr'{e) ■=^sgn{FF') d^r Ai{r^ iplg^^^py^p{r}(p[,g^^F'y]F'ir} , (A.l) 
the rates change to 



rr' ci 



X 

k'j 
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rr' (Ti 



X 

i'j 



rr' cTi 



where Ea—Ef, := Eat , and it was assumed that the tunnelhng and the reflection processes 
take mainly place close to the leads, which justifies to drop the position dependence of 
the electron operators. 

More details on this can be found in the appendix of [8], where it is also nicely explained 
how some considerations about the main contributions under the integrals in (lA.ip and 
the fact that the SWCNT is unpolarised with respect to the L- and -R-bands allow to 
set $/,.r'(e) = Srr'^i, and actually in analogue Rirr'{e) = Srr'Ri- Then it is rather easy 
to carry out the integration / de dt2- 

Any formulary tells us for some real function G{e): 

Re(^j deG(e) ^ dts e^^^^-^)*^^ = 7ihG{E) , 



lm( I deG{e) I dta e^^^^"^)*^ ) = ±n if de 



/ ^-oo 



e-E 



where ^ denotes a principal part integration 
So finally: 

TXL 



r (Ji 



X 

k'j 



D„.(E,,) ME,,) + 4 ( £y ^f^) 



r ai 



X 



A„(B..) (1 - ME,.)) - ai ( £de A.^U^-ZK^)) ^ 



(A.2) 
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